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Abstract. We have investigated the turbulent mean-field dynamo action in protoneutron stars that are subject to convective 
and neutron finger instabilities during the early evolutionary phase. While the first one develops mostly in the inner regions of 
the star, the second one is favored in the outer regions, where the Rossby number is much smaller and a mean-field dynamo 
action is more efficient. By solving the mean-field induction equation we have computed the critical spin period below which 
no dynamo action is possible and found it to be ~ Is for a wide range of stellar models and for both axisymmetric and 
non-axisymmetric magnetic fields. Because this critical period is substantially longer than the characteristic spin period of very 
young pulsars, we expect that a mean-field dynamo will be effective for most protoneutron stars. The saturation dipole field 
estimated by making use of the model of "global" quenching fits well the pulsar magnetic fields inferred from the spin-down 
data. Apart from the large scale magnetic field, our model predicts also a generation of small scale fields which are typically 
stronger than the poloidal field and can survive during the lifetime of pulsars. Extremely rapidly rotating protoneutron stars 
(P ~ 1 ms) may have the dipole field ~ (3 - 6) X 10 14 G. 
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1. Introduction 

The origin of pulsar magnetism is a subject of debate for 
decades. In a simple magnetic dipole braking model, the po- 
lar field strength inferred from observational data can reach 
~ 5 x 10 13 G. One possibility is that the magnetic field of 
a progenitor star is amplified by many orders of magnitude 
because of the conservation of the magnetic flux during the 
collapse stage (Ginzburg 1964, Woltier 1964). However, the 
"fossil field" hypothesis despite being seemingly attractive and 
plausible meets a number of problems. For instance, the pro- 
genitor star should possess a sufficiently strong magnetic field 
that does not agree with observational data (see discussion in 
Thompson & Duncan 1993 for more details). 

Another possibility is a turbulent dynamo action that can 
amplify the magnetic field during first ~ 30 — 40s of a neutron 
star life when the star is the subject of hydrodynamic instabil- 
ities. Hydrodynamic instabilities in protoneutron stars (PNSs) 
are driven by either lepton gradients which result in the so- 
called "neutron-finger instability" (Bruenn & Dineva 1996), or 
by negative entropy gradients which are common in simula- 
tions of supernovae (Bruenn & Mezzacappa 1994, Rampp & 
Janka 2000) and in models of PNSs (Keil & Janka 1995; Keil, 
Janka & Miiller 1996; Pons et al. 1999). The latter instability 
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is usually referred to as the "convective instability". The nature 
of instabilities in PNSs has been considered by a number of 
authors (Grossman, Narayan & Arnett 1993, Bruenn & Dineva 
1996, Miralles, Pons & Urpin 2000). Turbulent motions caused 
by instabilities in combination with rapid rotation that seems to 
be almost unavoidable in PNSs make turbulent dynamo one of 
the most plausible mechanism responsible for the pulsar mag- 
netism. 

Recent measurements of X-ray spectra of some pulsars 
have started to provide a closer look at the magnetic field at 
the neutron star surface. These observations indicate that the 
pulsar magnetic field may often have a fine structure near the 
surface apart from the global magnetic structure inferred from 
the spin-down data. For instance, the absorption features in the 
spectrum of IE 1207.4-5209 allow to estimate a surface mag- 
netic field as — 1.5 x 10 14 G (Sanwal et al. 2002), that is in con- 
trast with the dipolar magnetic field estimated from the spin- 
down rate of this pulsar (~ (2 — 4) x 10 12 G (Pavlov et al. 
2002)). Becker et al.(2003) have found an emission line in the 
X-ray spectrum of PSR B 182 1-24 that could be interpreted as 
cyclotron emission from the pulsar's polar cap. The line can be 
formed in a magnetic field ~ 3 x 10 n G, approximately two 
orders of magnitude stronger than the dipolar magnetic field. 
Kargaltsev, Pavlov & Romani (2004) have reported a marginal 
detection of the emission line in the spectrum of the millisec- 
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ond pulsar J0437-4715 that can be interpreted as a cyclotron 
line from the magnetic sport with the field ~ 7 x 10 8 G on 
the neutron star surface. Haberl et al.(2003) found evidence for 
a proton cyclotron line in the spectrum of the isolated pulsar 
RBS 1223. This line corresponds to the local magnetic field 
at the stellar surface of approximately one order of magnitude 
stronger than the dipole field. All these measurements provide 
evidence that the magnetic field likely has a complex structure 
in neutron stars with small scale fields being stronger than (or 
comparable to) the large scale field responsible for the secular 
spin-down of pulsars. 

Observations of radio emitting pulsars also support the 
idea that neutron stars apart from the large scale field (likely, 
dipoles) may have small scale magnetic structures at the sur- 
face. Recently, Gil & Mitra (2001) and Gil & Melikidze (2002) 
have argued that the formation of a vacuum gap in radio pul- 
sars is possible if the magnetic field lines near the polar cap 
have a small curvature <~ 10 5 cm and the field is very strong, 
B s ~ 10 13 G, irrespective of the magnetic field measured from 
the spin evolution. Furthermore, the presence of a strong mag- 
netic field with a small curvature can account for the radio 
emission of many radiopulsars that lie in the pulsar graveyard 
and should be radio silent (Gil & Mitra 2001). Analysing drift- 
ing subpulses observed in many pulsars, Gil & Sendyk (2003) 
found that their behavior is consistent with the vacuum gap 
maintained by a strong spot-like magnetic field. 

The growing number of evidences for a complex structure 
of the magnetic field suggests that this may represent a gen- 
eral phenomenon in neutron stars. The presence of a large scale 
field accompanied by small scale magnetic structures can nat- 
urally be understood if the field is generated by the turbulent 
mean-field dynamo. The duration of an unstable stage in PNSs 
(<~ 30 — 40 s) is sufficient for dynamo to reach a saturation level 
because the period of rotation of PNSs and the dynamo growth 
time are typically much shorter (Thompson & Duncan 1993). 
The magnetic fields generated by turbulent dynamo in convec- 
tive PNSs will be frozen in the crust that starts to form almost 
immediately after convection stops. Since the crustal conduc- 
tivity is high both the large and relatively small scale fields can 
survive in neutron stars during a long time comparable to the 
lifetime of pulsars (Urpin & Gil 2004). 

Recently, it has been shown by Bonanno, Rezzolla & Urpin 
(2003) that turbulence can drive both small and large scale dy- 
namos in PNSs. This result is in contrast to the previous con- 
clusion by Thompson & Duncan (1993) that only small scale 
dynamo can be operative in neutron stars. The reasoning of 
these authors was based on the assumption that the whole PNS 
is convectively unstable with turbulent velocity vt ~ 10 8 cm 
s _1 . An efficiency of the mean-field dynamo can be charac- 
terized by the Rossby number, Ro = P/tt, where P is the 
period of rotation and tt is the turnover time of turbulence. If 
P ~ 10 — 100 ms that is likely typical for young neutron stars 
(Narayan 1987) and tt ~ 1 ms that is typical for the convec- 
tively unstable region of PNSs then Ro ~ 10 — 100 and the 
influence of rotation on the turbulence is therefore weak. As a 
consequence, Duncan & Thompson (1992) and Thompson & 
Duncan (1993) concluded that the mean-field dynamo does not 
operate in PNSs except those rotating with the period <~ 1 ms. 



The model of pulsar magnetism qualitatively similar to that 
by Thompson & Duncan (1993) has been outlined by Wheeler 
et al. (2000) and Wheeler, Meier & Wilson (2002) who pointed 
out the importance of differential rotation that seems to be al- 
most unavoidable in PNSs. Their model suggests that a very 
strong toroidal magnetic field can be generated by differen- 
tial rotation in a newly born pulsar. These authors assume that 
a poloidal field of ~ 10 12 G could arise simply from flux- 
freezing if the precollapse progenitor core has a field strength 
comparable to that of a magnetized white dwarf, ~ 10 8 G. 
Mention that this point seems to be in contradiction with obser- 
vational data and has been criticized by Thompson & Duncan 
(1993). The poloidal field can be amplified further by differ- 
ential rotation to produce a strong toroidal field. The toroidal 
field is growing linearly with time, and the growth is limited by 
buoyancy which operates to expel the field from the site where 
it is generated. According to the authors, the toroidal field at the 
buoyancy limit can be as strong as ~ 10 16 G. Generally, such 
a strong field can be sufficient to suppress convection in proto- 
neutron stars (Miralles, Pons & Urpin 2002). Note, however, 
that the turbulent magnetic diffusion caused by convection can 
essentially decrease the saturation toroidal field as it will be 
seen from our calculations and, most likely, the generated field 
will not suppress convective motions. 

In fact, the picture is more complex since two essentially 
different instabilities may occur in PNSs (Bruenn & Dineva 
1996, Miralles, Pons & Urpin 2000). Convection is presum- 
ably connected to the entropy gradient, whereas the neutron- 
finger instability is more relevant to a negative lepton gradi- 
ent. The neutron-finger instability is the astrophysical analogy 
of the salt fingers that exist in terrestrial oceans. Physically, 
a fluid parcel perturbed downward in the PNS can thermally 
equilibrate more rapidly with the background but find itself 
lepton-poorer and denser and, therefore, subject to a downward 
force that would amplify perturbations. The neutron-finger in- 
stability is a sort of the so called "doubly diffusive" instabilities 
that occur due to dissipative processes like viscosity, thermal 
and lepton diffusivity. These processes are rather fast in PNSs 
and, therefore, the secular neutron-finger instability can grow 
rapidly. The estimated growth time in the neutron-finger unsta- 
ble region, ~ 30 — 100 ms, is only a couple orders of magni- 
tude longer than the growth time of extremely rapid convection 
(Miralles, Pons, & Urpin 2000), thus yielding a mean turbulent 
velocity — (1 - 3) x 10 6 cm/s. 

Typically, the convectively unstable region is surrounded 
by the neutron-finger unstable region, the latter involving there- 
fore a larger portion of the stellar material. As a result, not the 
whole PNS is convectively unstable in contrast to the assump- 
tion made by Thompson & Duncan (1993). Instabilities first 
develop in the outer layers containing ~ 30% of the stellar 
mass, The two unstable regions move towards the inner parts 
of the star and, after ~ 10 s, more than 90% of the stellar mass 
is hydrodynamically unstable. At this stage the stellar core has 
become convectively unstable but it is still surrounded by an 
extended neutron-finger unstable region. In <~ 20 s that follow, 
the temperature and lepton gradients are progressively reduced 
and the two unstable regions begin to shrink, leaving the outer 
regions of the star. After ~ 30 s, most of the PNS is stable and 
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the instabilities disappear completely after ~ 40 s (Miralles, 
Pons & Urpin, 2000). The Rossby number is actually large in 
the convectively unstable region, Ro ~ 10 — 100, and the mean- 
field dynamo likely does not work in this region. However, this 
is not the case in the neutron finger unstable region where the 
turbulent velocity is much slower than in the convectively un- 
stable region, and Ro ~ 1 (Bonanno, Rezzolla & Urpin, 2003). 
Therefore, turbulence can be strongly modified by rotation in 
the neutron finger unstable region, and this favors the efficiency 
of a mean-field dynamo. Of course, in both regions turbulent 
magnetic fields can also be generated by small scale dynamo. 

In the present paper, we consider the mean-field dynamo 
action in PNSs in more detail. We extend our study to the case 
of non-axisymmetric fields that are of particular interest in pul- 
sars. The main goal of this paper is to show that the mean- 
field dynamo can operate in a wide range of the parameters 
of PNSs and generate the field of the strength comparable to 
that of "standard" pulsars. The paper is organized as follows. 
In Section 2, we consider the basic equations governing the 
mean-field dynamo and discuss the properties of convection in 
PNSs. The numerical results are represented in Section 3, and 
our findings are summarized in Section 4. 

2. Basic equations 

To investigate the efficiency of a mean-field dynamo action, 
we model the PNS as a sphere of radius R with two substan- 
tially different turbulent zones separated at R c < R. The inner 
part (r < R c ) corresponds to the convectively unstable region, 
while the outer one (R c < r < R) to the neutron-finger un- 
stable region. The boundary between the two regions moves 
inward on a timescale comparable to the cooling timescale (i.e. 
<~ 1 — 10 s) that is much longer than the turnover time for 
the both unstable zones. The energy of turbulence is generally 
non-stationary as well, developing rapidly soon after the col- 
lapse, reaching a quasi-stationary regime after a few seconds, 
and then progressively disappearing. However, changes take 
place on the cooling time scale and, therefore, all parameters 
of turbulence can be treated in quasi-steady approximation. In 
this case, the mean-field induction equation for a turbulent PNS 
can be written as 

OB 

— = V x(vxB + aB) - Vx(i)VxB), (1) 

where rj is the turbulent magnetic diffusivity, a is a pseudo- 
scalar measuring the efficiency of the dynamo (the "a- 
parameter"), and v is the velocity of a mean fluid motion. 
Boundary conditions for the magnetic field need to be speci- 
fied at the stellar surface, where we impose vacuum boundary 
conditions, and at the center of the star, where we impose the 
vanishing of the toroidal magnetic field. 

We assume that differential rotation is the only large scale 
motion in PNSs, and v = fl(r) x r. From theoretical model- 
ing and simple analytic considerations it is commonly accepted 
that core collapse of a rotating progenitor leads to differential 
rotation of a newly born neutron star (Zwerger & Miiller 1997; 
Rampp, Miiller & Ruffert 1998; Dimmelmeier, Font & Miiller 
2002) mainly due to conservation of the angular momentum 
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during collapse. A recent study on evolutionary sequences of 
rotating PNSs (Villain et al. 2004) shows that the typical scale 
on which the angular velocity changes is in the range w 5 — 10 
km. We consider two possible models of differential rotation in 
PNSs, a shellular rotation with 

n(r) = nW+(l) 2 nS ) fc . (2) 

and cylindrical rotation with 

n oo = n S + (£) an 3. w 

where r and s are the spherical and cylindrical radii, respec- 
tively. Note that O(r) can generally depend on the time be- 
cause the PNS cools down and characteristics of turbulence 
change during the unstable stage. Turbulent transport of the 
angular momentum takes place on the diffusive timescale, ~ 
R 2 /v t tt, that is short compared to the cooling time of PNSs. 
Therefore, differential rotation is likely constrained to the prop- 
erties of turbulent motions at each instant of time, and f2(r) fol- 
lows adiabatically the PNS thermal and chemical evolution. We 
will neglect these relatively slow changes in our quasi-steady 
modeling. This can well be justified because the main conclu- 
sions of our study are qualitatively the same for any considered 
rotation law including rigid rotation. 

As it was mentioned, all of the PNS undergoes turbulent 
motions but with properties that are different in the inner parts 
(0 < r < R c ), where fast convection operates, from those 
in the outer parts (R c < r < R), where the neutron-finger 
instability operates. To model this in a simple way, we assume 
the relevant physical properties of the two regions to vary in 
a smooth way mostly across a thin layer of thickness AR = 
0.025i?. More precisely, we express rj as 

V = Vc + (Vnf - Vc) {1 + erf[(r - R C )/AR]} /2 , (4) 

where r\ c and r] n f are respectively the turbulent magnetic dif- 
fusivities in the convective and neutron-finger unstable zones, 
and erf is the "error function". Using equation (4), we have 
i] w t] c if R c — r 3> AR (convective zone), while 77 w rj n f 
if r — R c 3> AR (neutron finger unstable zone). In the both 
unstable zones, the turbulent magnetic diffusivity can be esti- 
mated as vt£t 1 3 where vt and It are the corresponding turbu- 
lent velocity and length-scale. The turbulent length-scales are 
likely comparable in the both zones and ^1 — 3 km. The turbu- 
lent velocities are, however, different with much larger velocity 
in the convective zone (see Bonanno, Rezzolla & Urpin 2003). 
Therefore, rj n f <C r\ c in PNSs, and we choose t] n j jr\ c — 0.1 in 
calculations. 

Similarly, we have modeled the a-parameter as being neg- 
ligibly small in the convectively unstable region and equal to 
a n f in the neutron-finger unstable region, i.e. 

a(r,6) = a nf cose {1 + erf[(r - R C )/AR}} /2 , (5) 

where the angular dependence chosen in this expression is 
the simplest guaranteeing antisymmetry across the equator. 
Usually, we assume a n f =const in calculations. However, the 
a n f -parameter is likely not constant in the neutron finger un- 
stable region because differential rotation or inhomogeneity of 
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Fig. 1. The sketch of the radial dependence of F(r) for the 
three models used in our calculations. 



turbulence can generally produce some non-uniformity in a. 
For instance, a non-uniformity can originate from overshooting 
caused by convection that smooths the boundary layer between 
the zones. To understand how much our conclusions are sensi- 
tive to the choice of a, we performed some calculations with a 
radially dependent a n f(r) — a n foF(r) where F(r) is plotted 
in Fig.l. 

We recall that in a rotating turbulence with length-scale 
It and moderate Rossby number, a„/o ~ Q£ T W ln(pv^) 
(Riidiger & Kitchatinov 1993). In PNSs, however, the pressure 
is mainly determined by number density of degenerate neu- 
trons and, therefore, the density length-scale is comparable to 
the pressure one L and to the maximal length-scale of the in- 
stabilities. This introduces a great simplification since we have 
then a n fo ~ QL. 

3. Numerical results 

The induction equation (1) with CI, -q and a given by equations 
(2)-(3), (4) and (5) has been solved with a numerical code em- 
ploying finite-difference techniques for the radial dependence 
and a polynomial expansion for the angular dependence. In par- 
ticular, we use the representation (Radler, 1973) 

B = B„ + B t , 



B p = V x A t 



(6) 



with 



B t = r x W, A t = -rx V$ (7) 
and the scalar functions ^ and $ have the following expansion 



* = E E < 

n—l m——n 



* = E E 



-i\„ 



-i\„ 



^ nm (r)Y™(6,< 



(8) 



n— 1 m- 



The induction equation thus decouple inamx n-dependent set 
of coupled equations for $ and $ which can be solved by trun- 
cating the order n of the harmonics for a given value of m (fur- 
ther details can be found in Radler, 1973, Bonanno et al, 2002). 
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Fig. 2. Critical period as a function of the differential rotation 
parameter q sp h for F(r) = 1. The two pairs of curves refer 
to different values of R c , with the solid parts corresponding 
to a stationary dynamo and the dashed parts to an oscillatory 
dynamo. The number m corresponds to different values of the 
azimuthal wavenumber in a polynomial expansion for the mag- 
netic field. 



Non axisymmetric modes, i.e. those with m ^ are waves 
traveling in azimuthal direction. The field configuration of a 
non-axisymmetric mode rotates like a rigid body with angular 
velocity \ mn /m. The simulations reported here use 30 spheri- 
cal harmonics and about 40 grid points in the radial direction in 
order to ensure convergence. Boundary conditions are chosen 
in order to guarantee the regularity at the center and vacuum 
outside r = R. We have solved the induction equation (1) to 
determine the critical value ao corresponding to the marginal 
stability of the dynamo. The seed magnetic field will grow if 
&nfo > ao and decay if a„/o < «o- Since a„/o ~ the 
critical value a effectively selects a critical value for the spin 
period, P = 2irL/a.o, such that magnetic field generation via 
a mean-field dynamo action will be possible only if the stellar 
spin period is shorter than the critical one. The different types 
of dynamo can be distinguished according to whether the gen- 
erated field exhibits periodic oscillations (oscillatory dynamo, 
dashed lines in Fig.2) or not (stationary dynamo, solid lines) 

In Fig.2, we plot the critical period Pq for the case of the 
shellular rotation (2). The critical period is shown as a function 
of the parameter q sph = 0^/(0^+0^) that characterizes 
differential rotation. Note that q > and q < correspond to 
situations in which the stellar surface rotates faster and slower 
than the center, respectively. Calculations have been done for 
rj n f = 10 11 cm 2 s _1 and F(r) = 1 (constant a-parameter in 
the neutron finger unstable region). The different types of dy- 
namo can be distinguished according to whether the generated 
field exhibits periodic oscillations (oscillatory dynamo, dashed 
lines in Fig.2) or not (stationary dynamo, solid lines). 

As shown in Fig.2, a stationary dynamo dominates the ax- 
isymmetric magnetic field generation process for \q\ < 1, while 
an oscillatory dynamo is more efficient for 1 < \q\. These two 
regimes correspond to a 2 -dynamo and afi-dynamo, respec- 
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Fig. 3. Critical period as a function of the differential rotation 
parameter q sp h for F(r) given by the model 2 in Fig.l. Other 
parameters are the same as in Fig. 2. 



Fig. 4. Critical period as a function of the differential rotation 
parameter q sp h for F(r) given by the model 3 in Fig.l. Other 
parameters are the same as in Figs. 2 and 3. 



tively. Given the large differential rotation required for a gen- 
eration of the axisymmetric oscillatory magnetic field, it may 
be difficult to achieve in practice. Hence, the a 2 -dynamo ap- 
pears to be the most likely source of axisymmetric magnetic 
field generation via dynamo processes in PNSs. The situation 
is, however, different for the generation of a non-axisymmetric 
magnetic field that is of interest for PNSs since the observed 
pulsars have non-axisymmetric fields. 

The critical spin period found in our calculations is in gen- 
eral rather long. For instance, a mean-field dynamo will de- 
velop if P < 1 s when R c /R = 0.6 and if P < 1.4 s when 
R c /R — 0.3 if a PNS rotates rigidly (i.e. q sph = 0). If the 
PNS rotates differentially, the critical period is typically re- 
duced. The difference, however, is not large except the case of 
the non-axisymmetric field (m = 1) in a star with \q sp h\ > 1 
when the difference can reach ~ 40%. As a result, only PNSs 
with P > 1 s will not be subject to a turbulent mean-field dy- 
namo action. Such slow rotation rates should be rather difficult 
to achieve if angular momentum is conserved during the col- 
lapse to a PNS. We expect, therefore, that a turbulent mean- 
field dynamo will be effective during the initial stages of the 
life of most PNSs. 

In Fig. (3) and (4), we plot the same as in Fig. (2) but for 
F(r) represented by the lines 2 and 3 in Fig.(l), respectively. 
The behavior of Pq is qualitatively unchanged and the only dif- 
ference is that there is no stationary dynamo regime for a gen- 
eration of the axisymmetric field in the case R c /R = 0.6 if 
the function F(r) is decreasing to the surface (the model 3 in 
Fig.(l)). Like the case F(r) — 1, the critical period is longer 
for R c /R = 0.3 compared to R c /R — 0.6. The difference 
is due to the fact that a PNS with a more extended neutron- 
finger unstable region can rotate proportionally more slowly 
while maintaining the same dynamo action. The values of Po 
are typically by a factor ~ 2 smaller than for a dynamo action 
with F(r) = 1 implying that the onset of dynamo is mainly 
characterized by the average value of a in the unstable region 
rather than a particular dependence of a on r. 
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Fig. 5. Critical period as a function of the differential rotation 
parameter q cy i for F(r) given by the models (1) and (2) in 
Fig. 1. for R c /R = 0.6. 



The above discussion for the shellular model Eq.(2), re- 
mains qualitatively similar for a cylindrical rotational law 
Eq.(3) for q cy i = + £2^). Since in this case the 

differential rotation is not important the basic dynamo action 
is that of an a 2 -model. If instead q cy i > 1, the dynamo mech- 
anism would produce an oscillating field of an af2-type. The 
general dependence of the critical period on the differential ro- 
tation parameter q cy i is shown in Fig.(|5j for R c = 0.6R. Note 
that, for small q cy i, the m = 1 modes have a greater critical 
periods than the m = modes, and thus the critical value of 
ct n fo for non-axisymmetric modes is smaller than for the ax- 
isymmetric case. As a consequence, non-axisymmetric field is 
more easily excited than the axisymmetric one in this case. 

In Fig.l|6j, we show the toroidal (B^) and poloidal (B p ) 
magnetic fields of a typical PNS model for the cylindrical ro- 
tation law. Note that the both components are generated in 
the outer neutron-finger unstable region but turbulent diffusion 
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Fig. 6. Toroidal (B$) and the poloidal {B p ) magnetic field lines 
for different pulsational times for m = and q w 2 in Fig.Q. 
Solid and dashed contours correspond to positive and negative 
values, respectively. The dot-dashed line marks the position of 
R c = 0.6R. 



transports the magnetic field also in the inner convectively un- 
stable region. This field is, however, considerably weaker. Our 
calculations show that if \q\ < 1 and the field is generated by 
the a 2 -dynamo then B^/B p ~ 10, while B^/B p ~ 100 - 200 
if I q I > 1 and the ail-dynamo generates the magnetic field. 
Both results suggest that the internal magnetic fields in neu- 
tron stars could be substantially stronger than the observable 
surface fields. The toroidal magnetic field tends to concentrate 
near the polar regions whereas the poloidal one is more evenly 
distributed in latitude. Note that, in the case of an a 2 -dynamo, 
the generated field propagates more efficiently into the convec- 
tive zone (see Bonanno, Rezzolla and Urpin 2003), 

4. Discussion 

The PNS is subject to two substantially different instabilities, 
with a convective instability operating in the inner region of 
the star and a neutron-finger instability being more efficient 
in the outer region. The turbulent motions are more rapid in 
the convective zone, where the Rossby number is large but the 
a-parameter, characterizing the mean-field dynamo action, is 
likely small. In the neutron-finger unstable region, on the other 
hand, the turbulent turnover time is considerably longer, the 
Rossby number small, and the a-parameter can be sufficiently 
large to drive a mean-field dynamo. 

The a-parameter depends the stellar rotation being larger 
for rapidly rotating stars. Our simulations show that even rel- 
atively slowly rotating PNSs can be subject to a dynamo ac- 
tion, with the a 2 -dynamo being the most efficient mechanism 



of generation for both axisymmetric and non-axisymmetric 
fields if differential rotation is not extremely strong. The cal- 
culated critical value of the spin period that determines the 
onset of dynamo in PNSs is P ~ 1 s for a wide range of 
models. This value is essentially larger even than the charac- 
teristic spin period of very young pulsars as inferred from ob- 
servations (~ 50 - 100 ms, Narayan 1987) but likely PNSs 
can rotate even faster. As a result, a turbulent mean-field dy- 
namo can be effective in the early stages of the life of most 
PNSs. The critical period is not very different for axisymmet- 
ric and non-axisymmetric fields and likely both these magnetic 
configurations can be generated by the PNS dynamo. The gen- 
eration of a non-axisymmetric component is the key point for 
the magnetism of pulsars since they have a substantial non- 
axisymmetric field. It must be stressed that, as the neutron- 
finger unstable region shrinks towards the surface, the differ- 
ence between axisymmetric and non-axisymmetric critical pe- 
riods, disappears as it can be noticed, for instance, in Fig. 2 
and Fig. 3 thus providing an opportunity to generate a non- 
axysimmetric field in a very general framework. 

The PNS dynamo can operate both in oscillatory and sta- 
tionary regime. A non-axisymmetric field is always a wave 
traveling in the azimuthal direction. On the other hand, axisym- 
metric field can develop in both the oscillatory and stationary 
regimes, depending on the profile of a in the neutron finger un- 
stable region and the thickness of this region. The period of dy- 
namo oscillations can be estimated as ~ Q.1R 2 /rj n f ~ R 2 /r] c 
and, in our model, is ~ Is. Note that, for the pulsar mag- 
netism, there is no particular difference in which regime, oscil- 
latory or stationary, the field is generated because a formation 
of the crust starts almost immediately after the instability stops, 
and the generated field should be frozen-in a highly conductive 
crustal matter. 

The unstable stage lasts ~ 30 — 40 s in PNSs, and likely this 
time is sufficient for the dynamo to reach a saturation level. We 
estimate the field strength in saturation by making use of the 
model of a "global" quenching. In accordance with this model, 
a generation of the magnetic field decreases the a-parameter, 
and the simplest possible expression can be proposed for the 
non-linear a-parameter in the neutron finger unstable zone, 



a n i{B) 



a n fo 



B 2 /Bl q 



(9) 



where B is the characteristic value of the generated field, and 
B eq is the equipartition magnetic field determined by equating 
the kinetic and magnetic energy of turbulence. Then, the satu- 
ration magnetic field, B sat , can be estimated if we assume that 
otni{B S at) is equal to the critical value ao that determines the 
marginal dynamo stability. We have from this condition 



B„ 



B e 



P 



1, 



(10) 



If Pq > P > Po/2, the generated mean field is weaker than the 
turbulent magnetic field, B eq . On the contrary, the saturation 
field is stronger than B eq if P < Pq/ '2, and B sat can reach a 
rather large value if the PNS rotates very rapidly. Likely, how- 
ever, that the rotation period cannot be shorter than ~ 1ms 
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and, hence, the maximum field generated by dynamo in PNSs 
is around ~ 30B eq . In our calculations, the generated toroidal 
field is typically stronger than the poloidal field. Therefore, 
quenching is determined mostly by the strength of the toroidal 
field, and equation (10) provides an estimate of the saturation 
toroidal field. Since, in the case of an a 2 -dynamo, the poloidal 
field is approximately 5-10 times weaker, we obtain for the sat- 
uration poloidal field 



S paot w(0.1-0.2)Be,y^-l. (11) 

The poloidal field is typically weaker than (or comparable to) 
the equipartition field. The only exception is very rapidly ro- 
tating PNSs with the period ~ 1 — 3 ms where the poloidal 
field can be a factor few stronger than small-scale turbulent 
fields. In the majority of pulsars, however, we can expect that 
a large scale poloidal (for example, dipolar) field should not 
be stronger than the small-scale fields if their magnetism is 
caused by the mean-field dynamo action. Note that because of 
a high conductivity, turbulent magnetic fields with the length- 
scale > 1 km can survive in the neutron star crust for a very 
long time comparable to the lifetime of pulsars, <~ 100 Myr 
(Urpin & Gil 2004). 

Turbulence is non-stationary in the both unstable zones 
of PNSs, and therefore B eq in equation (9)-(ll) varies with 
time. It rises very rapidly soon after collapse, reaches some 
quasi-steady regime, and then goes down when the tempera- 
ture and lepton gradients are smoothed (after ~ 30 — 40 s). The 
timescale required for fluid to make one turn in a turbulent cell 
with the length-scale It can be estimated as tt ~ it£t/vt- 
This timescale varies with time as well, but is typically much 
shorter than the characteristic cooling timescale of the PNS, 
t coo i, except the very late phase when gradients are smoothed 
and instabilities are less efficient. Therefore, turbulence can 
be treated in a quasi-steady approximation during almost the 
whole convective evolution of the PNS except the late stage. 
We can estimate B eq during the quasi-steady regime as ~ 10 16 
G in the convective zone, and <~ (1 — 3) x 10 14 G in the neu- 
tron finger unstable zone (see Urpin & Gil 2004). However, 
the temperature and lepton number gradients are progressively 
reduced as the PNS star cools down and, therefore, the tur- 
bulent velocity decreases as well. As a result, the strength of 
small-scale magnetic fields generated by turbulence also de- 
creases compared to the maximum value, but the turnover time 
of turbulence increases. The mean-field dynamo as well as the 
small-scale one are still operative until the quasi-steady con- 
dition T coo i ^> tt, is fulfilled. We assume that this condition 
breaks down at some instant of time when tt becomes compa- 
rable to the cooling timescale: tt ~ t coo i. Then, the turbulent 
velocity at this instant is given by 

vt ~ • (12) 

Tcooi 

We assume that the final strength of the magnetic field gener- 
ated by both the mean-field and small-scale dynamo is deter- 
mined by B eq at the instant of time when the quasi-steady con- 



dition breaks down. Therefore, we have for the final strength of 
the equipartition magnetic field 

r- — TTy/Mpe T 

B eq ~ V ■iirpvT ~ • (13) 

^cool 

The final strength of the generated small-scale field turns out to 
be the same for both unstable zones. For the largest turbulent 
scale, i T = L ~ 1 - 3 km, estimate (13) yields B eq ~ 3 x 10 13 
G if Tcooi is of the order of a few seconds. Using this estimate of 
B eq , we can conclude that the strength of a large-scale poloidal 
field generated by the mean-field dynamo (equation (8)) is in a 
good agreement with the observed magnetic fields of the ma- 
jority of pulsars. For example, the generated poloidal field is 
— (1 - 2) x 10 13 G if the star rotates with the period — 100 
ms. Note that extremely rapidly rotating PNSs (P <~ 1 ms) 
may possess a very strong poloidal field <~ (3 — 6) x 10 14 G 
comparable to that of magnetars. 
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